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*rH , Abstract 

a: 

o 1 We investigate the SM Higgs decays into heavy quarkonia J/ip and 

T in both color-singlet and color-octet mechanisms. It is found that 
in J/yb production the contribution through color-octet processes over- 
.£h ! whelm that in color-singlet processes all the intermediate region Mz < 

Mho < 2M W and the fraction ratio is comparatively large, but in T 
production the contribution through color-octet mechanism is negligible. 
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Despite the success of the Standard Model (SM) |l| in describing elementary particle 
phenomena at electroweak energy scale and below, the theory has some intrinsic shortcomings 
which may not be overcome in the context of the theory itself. For example, the hierarchy 
problem existed in the gauge symmetry breaking. The mechanism of spontaneous symmetry 
breaking introduced in SM to allow some of the gauge bosons and fermions to acquire mass 
while keeping the theory renormalizable leads to a neutral scaler particle, the Higgs boson 
H°. While a worldwide experimental effort has been made, this Higgs particle has not yet 
been discovered. 

Though the SM predicts the existence of the H° with unspecified mass, the couplings of 
it with other fundamental particles are well defined. Consequently the H° production cross 
section in various channels and the decay branching ratios of it are precisely predicted as 
a function of the Higgs boson mass. According to the LEP I data, based on the Bjorken 
process e + e~ — > H°Z*, the Higgs mass, M H o, should be heavier than 77.1 GeV 0. 

The hunt of the Higgs boson meets two problems, i.e., how to produce it and then how to 
detect it. Therefore, experimental researches at LEP include both H° production and decays, 
meanwhile the search range will expand with the operation of LEP at higher energy. After 
the operation of the Large Hadron Collider (LHC), almost all the possible energy region 
where Higgs may exist will be searched down. In this paper, we proceed out discussion 
in the intermediate energy range Mz < M^o < 2M W , which is the most difficult range in 
experiment to detect the Higgs boson. 

In theory, the study for the nature of the Higgs boson have been in highlight in particle 
physics ever since the establishment of SM. The probabilities of finding the H° and other 
Higgs bosons appearing in theories beyond the minimal SM in various collider facilities are 
evaluated || [Q . Among others of Higgs boson production processes, the Yukawa process, 
i.e., the bremsstrahlung production of the neutral Higgs boson H° from a heavy fermion, 
of the SM Higgs boson is an important channel to discover and to investigate this particle. 
In experiment, the analysises of the Yukawa process have led to some constraints on the 
theories in respect to the Higgs sector |j| . Among various decay processes of the Higgs boson, 
Fermionic decay channels of it are important channels. The observation of such decay modes 
may provide important information on fermion Yukawa couplings which associate with the 
symmetry breaking mechanism in the SM. The subprocess of the fermionic decay channels, 
that is quarkonium production in the Higgs boson decays that may be called the Extended 
Wilczek Processes (EWP), which helps to further identify the Higgs boson signals also gains 
much attention though there may be large theoretical uncertainties both due to the Quantum 
Chromodynamics (QCD) and relativistic corrections 0. Considering the fact that the low 
limit of the SM Higgs mass is up to about the mass of Z° boson's and the recent developments 
in quarkonium production theory, in this paper we reevaluate the EWP within the scope of 
Standard Model. 

In the past few years, two approaches in describing the heavy quarkonium produc- 
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tion, the fragmentation mechanism Q and the Non-Relativistic Quantum Chromodynamics 
(NRQCD) factorization ||, have greatly improved the theoretical studies in this area. The 
introduction of color-octet mechanism asks ones to reevaluate almost all the processes associ- 
ated with quarkonium production, in order either to test the new production mechanism, to 
further understand QCD or to assess the quarkonium production rates in a specific process. 
Of course, the EWP is also in the case. 

For charmonium and bottomonium production in the Higgs boson decay, the important 
color-singlet channels include H° — > V'j (electromagnetism), H° — > QQV (quark fragmenta- 
tion), H° — > bbV (photon fragmentation), and H° —>■ bbVgg (gluon fragmentation) processes 
as shown in Fig.l, and the process H° — > Vgg vanishes. Here, V represents a vector meson 
consisting of a QQ pair (Q is heavy quark b or c). 

The decay width of the process of Fig. 1(a) can be calculated straightforward making use 
of the standard formalism ||. 

r = ^Mi-^)<o|or( 3 Si)|o>, (i) 

where is the charge fraction of heavy quark, e c = | and e& = — ~; v is the characteristic 
quantity which sets the scale of spontaneous symmetry breaking in SU(2)i x U(l) theory 
and v ~ 246 GeV; £ y is defined as £ v = jf*-, M v and M H o are the masses of the bound 
systems and the SM Higgs boson; < 0|(9] / ( 3 S'i)|0 > is the matrix element of color-singlet 
four fermion operator of quarkonium, which is associated with the nonperturbative part 
of the quarkonium production and its value may be determined by quarkonium electronic 
decay widths or from the potential model. Because the Higgs H° couples preferentially to 
the heavy fermion, when M H o < 2M W , M w is the mass of W boson, the dominant decay 
channel in Higgs decays is H° — > bb even considering the virtual Z* and W* process [fTOfl . 
So, here after, we will constrain the H° coupling only to the quarks (b, c) in our discussion. 



r M = |4^o(i-V) 3 / 2 , ( 2 ) 

where rj = mb/M^o. From Eq. (|Tj) and Eq.(H), the fraction ratio 

_ r g o^j/^ 7 32TO Q/^(l n\ n J/j>(3g Mn ^ (o\ 

T H0 ^ b - b 729M H om^ (1 - Ar] 2 y/ 2 



_ T H o^ Tj 8na gr(l--. gj n|f?Tf 3^i n , 

= " 729MW(l-V)3/ 2 < ( > 

"72^^ 1 -*)- 1A < ' O fWl°>' (4) 

can then be easily obtained. The decay H° V'j with the subsequent decay V — > e + e~ 
from experimental point of view presents the clearest signal in searching the Higgs boson. 
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Unfortunately, the ratios of these decay rates to the decay rate H° — > 66 is too small (see 
below and Ref.0). 

Because the minimal value of the mass of H° is over 77.1 GeV, the decay widths, therefore 
the fraction ratios of the process shown in fig. 1(b) can be evaluated using the universal quark 



fragmentation functions [12 



_ T(H^ccJ/jj) m\ 

Rji* = T{H o - b - b) ~ 2 x p ~m x ( 5 ) 



where P c ^j/^ is the universal fragmentation probability, and 



16 /n , 2 < 0|Ci^( 3 5i)|0 > .1189 
243 30 



As Eq. (|S|) , (Q) , the fraction ratio of r H o_^x over ^ H°~>bl can be readily obtained, 

_ r(#° -> gr) 

= YiiP^W ~ 2 x P ^ T - (7) 

Here Pb^r, the fragmentation probability of b — > T, is in the same form as Eq.(^), 

16 /n , 2 < 0|Oj r ( 3 S'i)|0>.1189 c „ o , 
P fe _ x = ^a.(W ^ (-35- - 57in2). (8) 

The process as shown in Fig. 1(c) can be calculated using the standard method in com- 
puting the quarkonium production and decays and the result is simplified into the following 
form after neglecting the mass of the 6(6) quark (antiquark), except for in the Higgs-quark 
coupling vertex. 

dT{H° -> Vbb) _ a 2 m 2 b e 4 Q M H o < 0\OY( 3 Si)\0 > 



dz 1 dz 2 54M^v 2 7r z\z\ 



-[{z x z 2 - e v )(4 + 4) + 2z lZ2 ( Zl - l)(z 2 - 1)], (9) 



where z\ = 1 — fff 2 , z 2 = 1 — jjP- and their integration regions are r < z± < 1, r/z± < 
z 2 < 1 +r — z\. k, p, pi, p 2 are the momenta of the Higgs boson, bound state, and the final 
quarks appearing in Fig. 1(c). Because Eq.(|9|) is a differential equation, it can be used to 
evaluated the energy distribution of quarkonia, which is different from that of the two-body 
decay mode. After integrating over z\ and z 2 , the total decay width can be obtained and 
hence the fraction ratios. 

The decay width of the process H° — ► bbVgg as shown in Fig. 1(d) can be expressed as 

r(H°^bbg*,g*(ix)^ggV)= J d^T(H G -> bbg*(p,))P g ^ ggV . (10) 

M 2 

Here, P g *^ gg y is the decay distribution of the virtual gluon with virtuality /i, and it is defined 

as 

W - ggV) = 7r^P g ^ ggV . (11) 
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The expression of T(g* — ► ggV) can be found in Ref . |13| , the Pg*^ gg v, therefore, may be 
induced. 



l+r 



H 2 P(g* ->Vgg) = C v r J dx v J dx 1 f(x v ,xi,r), (12) 

where r = My/fi 2 , and 



lQgg < OIO^SQIO > lOag <0|O 1 T ( 3 5 1 )|0> 

^ 243vr 2 M] /v , ' 243vr 2 Mf ' 1 j 

The function / in Eq.flT2]) is of the form 

,/ \ (2 + x 2 )x 2 (2 + x 1 )x 1 

f{Xv,Xi;r) = — yry- + 



+ 



(2 - x JH ,)\l - Xl - r) 2 (2 - xj/^fil -x 2 - r) 2 

(xj/^-r) 2 -l 1 / 6(l + r-xj/^) 2 



'1 — x 2 — r) 2 (l — xi — r) 2 (2 — Xj/^) 2 ^ (1 — X2 — r) 2 (l — x\ — r) 2 



+ 2(1-^1(1-0 + 1 } 
(1 — X2 — r)(l — Xi — rjr r y 

where x, = 2^ / /i with i = V, g%, g 2 are the energy fractions carried by the quarkonia and 
two gluons in the g* rest frame, and then x 2 = 2 — xi —xj/^. The limits of the X\ integration 
in Eq. (Q) are 

^± = ^(2-x y ±V / 4-4r). (15) 

T(H° — ► bbg*(n)) is the same as Eq.(|S|) with just My being changed to Therefore, the 
fraction ratio H° — > bbVgg to if — > 66 then can be calculated using Eqs.flU) and (|10|). 

The heavy quarkonium production color-octet mechanism provides the process H° — > 
QQ( 3 Si,8)g — > Vg with soft hadrons, as shown in Fig.2(a), existing as the leading order 
process in a s . Although this process makes a less significant contribution to the J/ip produc- 
tion in Z° decays [|14] , comparing with other higher order processes, it may play a relatively 



more important role in the case of H° decays. On the other hand, after summing over the 
soft hadrons that accompany the quarkonium V and over the collinear hadrons that make 
up the gluon jet, the decay rate can be factored into the product of a two-body decay rate 
and a nonperturbative matrix element, and the two-body decay mode has a unique distri- 
butions in the kinematic variables of final states relative to that of three-body's. Another 
color-octet process which might play an important role in charmonium production in H° 
decays is shown in Fig. 2(b). 

Although the nonperturbative color-octet matrix elements in these color-octet processes 
are suppressed by v 4 with respect to that of color-singlet processes according to NRQCD 
velocity scaling rules, the strong coupling constant a s greatly enhanced them relative to the 
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corresponding electromagnetic processes. The calculation of the decay widths of these color- 
octet channels is similar to the corresponding color-singlet processes, with only the change 
of coupling constant, the color factor, and the nonperturbative matrix elements. Therefore, 
H° decay widths of the processes as shown in Fig.2 can be inferred from Eq.(|I|) and Eq.(H). 

F(H° -> VgX) = 0£v(l - £v) < 0|O 8 V ( 3 Si)|0 > . (16) 

Here, the X stands for the soft hadrons produced in the neutralization process of the color- 
octet (QQ)s evolving into the heavy quarkonium. The differential decay width of the Fig.2 (b) 
process reads as 

dT(H° -> VbbX) _ a 2 s m 2 b M H o < 0|£>^( 3 Si)|0 > 2 2 2 



dz x dz 2 12M 3 v 2 n z\z\ 



\{zxz 2 - H){z\ + z 2 2 ) + 2z 1 z 2 (z 1 ~ 1)0*2 - 1)],(17) 



where the integration variables z\ and z 2 change in the same regions as that in Eq.([§). In 
the above < 0|Q^( 3 Si)|0 > is the matrix element of four fermion color-octet operator. Its 
value may be determined by the fitting of theoretical predictions to the experimental data 
or from the lattice QCD calculations. Note that in the process of Fig.2 (a) the color-octet 
QQ^Sf>) and QQ( 3 Pj 8) ) have zero contributions to the H° decay widths. From Eqs.(||), 
(|16|), and fllTD the fraction ratios of quarkonium production over b-quark pair production 
then can be obtained. 



In the numerical calculation, we take 12 
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a s {2m c ) = 0.26, a s {2m b ) = 0.19, 
Mj/^p = 2m c , M r = 2m b , m b = 4.9 GeV, m c = 1.5 GeV, 

< OlOf^fSOlO >= 1.2 GeV 3 , < OlOj^Sx)^ >= 9.3 GeV 3 



< 0|Og / ^( 3 S'i)|0 >= 6.6 x 10~ 3 GeV 3 , < 0|O 8 T ( 3 Si)|0 >= 5.9 x 10~ 3 GeV 3 . (18) 

With the mass of SM Higgs boson changing from 65 GeV to 155 GeV, the relative weights 
and the magnitudes of the total fraction ratios of Rj/^p and R~c are shown in Fig. 3 and 
Fig. 4. There are multiple components to each of the color-octet and -singlet contributions 
in these figures. For the convenience in comparing the relative importance of those diagrams 
in Fig.l and Fig.2, we plot each component of the octet and singlet in Fig. 3 and Fig.4. 
Each figure of them in fact contains nine curves corresponding to six subprocesses of Fig.l 
and Fig.2, the sum of both octet and singlet subprocess, and the total. However, because 
the diversities of each component are so large that some curves overlap, which means the 
dominant contribution is almost the same as total. 

From the results shown in Fig. 3 and Fig.4, it is interesting to see that the heavy quarko- 
nium production in MS Higgs decays has two distinct features. First, generally, within the 
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color-singlet model the T production rate is smaller than that of J/^'s at collider facilities, 
but in the SM Higgs decays where they have a comparable fraction ratios because of the 
fact that the Higgs tends to couple with heavy fermions. From Fig. 3 and Fig.4 we can see 
that the process Fig. 1(b) always takes the leading contribution among all color-singlet ones 
no matter in J / if) production or T production. 

Second, the color-octet mechanism makes a dominant contribution in the J /if} production, 
however it makes a negligible contribution in the T production. This is because that for 
the intermediate mass Higgs the color-octet decay mode in Fig. 2(b) always overwhelms that 
in Fig.2(a), but in Fig. 2(b) gluon propagator in T production suppresses the decay width 
relative to the J / if) production. This unique character provides the EWP of J / if) production 
more reachable in experiment than T production if the color-octet mechanism does really 
work. In other words, if we would have measured a larger production rate in T production 
than in J /if) in the Higgs decays, we may conclude that the color-singlet mechanism is the 
dominant one in quarkonium production, otherwise, the other production mechanism, e.g., 
the color-octet one, must have taken part in. 

For the LHC luminosity of JC = 10 B p6 _1 there will yield 0(0.2-1 x 10 7 ) SM Higgs events 
in the mass range 155 GeV > M#o > 65 GeV with cm. energy y^s = 14 TeV ||16|| . From 
this the quarkonium events produced at the proton collider LHC can be readily estimated. 
For instance, at a mass of M H o = 100 GeV, from Fig. 3 and Fig.4 we can read that the Rj/^ 
and Rr to be about 10~ 4 and 2 x 10~ 5 , from Ref.[l6[] we get to know that the H — > bb 
branching ratio is about 80%. So, at this point the cross section times the branching ratio 
for pp — > H(— > J /if), T) + X will yield about 200 J/if> and 80 T events. 

In conclusion, after taking the color-octet mechanism into consideration the quarkonium, 
especially the charmonium, decay mode of SM Higgs may be detectable in the next generation 
collider, the LHC. That means the concerned quarkonium production process discussed in 
this paper may stand as an implement to the two-photon decay mode in identifying the SM 
Higgs in the intermediate mass region. 
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Figure Captions 



Fig.l. Color-singlet processes of H° decays into quarkonia, (a) electromagnetism process (b) 
quark fragmentation process (c) the photon fragmentation process (d) gluon fragmentation 
process. 

Fig. 2. Color-octet processes of H° decays into quarkonia, (a) the lowest order QCD process 
(b) gluon fragmentation process. 

Fig. 3. The total fraction ratios, relative weights of the different J ftp production mechanisms, 
and different subprocesses in H° decays. The dotted curve and the dashed curve illustrate 
the total fraction ratios of the color-singlet and color-octet mechanisms, respectively, the 
solid curve is a sum of them, and the subprocesses contributions are labeled a-d and A-B 
corresponding to the ones in Fig.l and Fig. 2. 

Fig. 4. The total fraction ratios, relative weights of the different T production mechanisms, 
and different subprocesses in H° decays. The curves in this figure are labeled the same as 
those in Fig.3. 
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